Lattice Approximation for Stochastic 
Reaction Diffusion Equations with One- Sided 

Lipschitz Condition 

Martin Sauer*^ and Wilhelm Stannat** 



^T) We consider strong convergence of the finite differences approximation in space for 

stochastic reaction diffusion equations with multiplicative noise under a one-sided 
Lipschitz condition only. We derive convergence with an implicit rate depending 
on the regularity of the exact solution. This can be made explicit if the variational 
solution has more than its canonical spatial regularity. As an application, spatially 
^ extended FitzHugh-Nagumo systems with noise are considered. 
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1 Introduction 



Oh 

Stochastic partial differential equations arise as a model in many scientific fields and the for 
applications much needed theory of numerical solutions to such equations is growing successfully. 
However, it is often the case that particular examples such as stochastic reaction diffusion 
equations, even with very simple structure of the nonlinear term, fail to fit in the assumptions 
of the majority of publications in this research area. In this article, we study such an SPDE of 
reaction diffusion type where the nonlinear drift term satisfies a one-sided Lipschitz condition 
only and is of the form 

v(t, = 6%v(t, + 0i (£, v(t, 0, w(t, 0) + b(£, v(t, 0,w(t, 0)v(t, £) 1 

for t > 0, £ G (0, 1), where rj(t,^) is some Gaussian white noise. This system consists of two 

variables v and w, where the former satisfies some semi-linear stochastic evolution equation 

and the latter is a kind of control variable coupled to v in an equation without diffusion and 
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no further independent external noise. The motivation for such a system comes from the field 
of neurobiology, in particular we consider the spatially extended stochastic FitzHugh-Nagumo 
system, modeling the propagation of the action potential in the axon of a neuron, see e. g. |2j [12] 
and Section 6 for details. 



The goal of this article is to prove /^-convergence including explicit error estimates of an 
easy to implement spatial approximation scheme. We consider the well-known finite difference 
method which has been studied by many authors, but to the best of our knowledge there is 
no convergence result which gives explicit rates if the nonlinear drift term is not Lipschitz 
continuous. Crucial references in the literature studying similar equations to ([I]) are Gyongy 
[5], where uniform convergence in probability without a rate is proven, as well as Pettersson & 
Signahl |15] for a similar result with multiplicative noise. Both results are obtained without 
a global Lipschitz condition. Moreover we want to mention a series of articles by Gyongy &; 
Millet |6| [7] , where the authors work in the variational approach and incorporate various spatial 
approximation schemes for SPDEs driven by a finite dimensional noise. However, equation 
([TJ is beyond their results since it excludes polynomially growing nonlinear ities. On the other 
hand, it is worth mentioning that there are existing results by Jentzen and Liu |13j which 
incorporate our one-sided Lipschitz assumption, but use a different approach, the spectral 
Galerkin approximation, to show path-wise convergence and ^-convergence with explicit rates, 
respectively. 



The SPDE ([I]) is studied within the variational approach, see e.g. |16| . and we only use the a 
priori known regularity of the exact solution to prove ^-convergence with an implicitly given 
rate in Theorem |3.1| This rate is given in terms of the regularity of the exact solution. If there 
is a better a priori information on the exact solution, this rate can be made explicit in terms of 
the approximation parameter. The proof is essentially based on uniform exponential a priori 
estimates for the approximating solutions obtained in Proposition 4.3 and the ltd formula for 
the variational solution. 



Most of the above mentioned articles concern fully discrete schemes, i. e. a discretization of 
the space and time domain. We only consider the spatial approximation part, which reduces 
the infinite dimensional problem to an SDE on a finite dimensional space, and then has to 
be solved with existing theory and the errors accumulate. Nevertheless, let us make a brief 
comment on the time discretization. One of the most popular approaches is the Euler-Maruyama 
scheme because of its simplicity and low computational effort. It is known that this scheme 
converges to the exact solution of the SDE for equations with globally Lipschitz continuous 
drift and diffusion coefficients. Of course, also our approximated equation fails to satisfy such a 
condition. For such equations with super-linearly growing drift and/or diffusion it was shown 
recently by Hutzenthaler, Jentzen &; Kloeden [9| that the Euler-Maruyama approximations 
diverge in a strong and weak sense. An obvious remedy for this problem is to use e. g. the 
implicit Euler scheme, which is known to converge, however it requires more computational 
effort than the explicit scheme. Again Hutzenthaler, Jentzen & Kloeden |10] suggest a so-called 
tamed Euler scheme that is explicit and computationally less expensive. It is also observed that 
the runtime as a function of the dimension grows linearly and quadratic for the tamed Euler 
and implicit Euler scheme, respectively. Of course, this favors the tamed Euler scheme for our 
high dimensional approximating problem. Thus, we propose solving ([I]) with a combination of 
finite differences and the tamed Euler scheme. 
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The article is structured as follows. In the next section, we describe the precise setting and 
assumptions on the coefficients of ([I]) and state the existence and uniqueness result for its 
corresponding abstract stochastic evolution equation in the sense of variational solutions. In 
Section 3 we introduce the approximation scheme and state the main result in Theorem |3.1| 
Sections 4 and 5 contain the proof of this theorem. The last part is devoted to the application of 
our results to a spatially extended FitzHugh-Nagumo system with noise. This model has been 
studied by Tuckwell [18J concerning the impact of noise on the generation of action potentials 
and the reliability of faithful signal transmission. Theorem 3.1 now allows to explicitly quantify 
the approximation error of the finite difference approximation. As an illustration for the use of 
our result, we state an estimator for the probability of propagation failure based on numerical 
observations as well as confidence intervals for this estimators. These depend on the statistical 
Monte-Carlo and the numerical approximation error and are explicitly constructed. 



2 Mathematical Setting and Assumptions 



Let (H, \\-\\h) = L 2 (0, 1) with scalar product denoted by (•, •) and consider the Laplacian as a 
linear operator on H, A : D{A) C H — > H, i. e. 

{Av)(0 :=^«(0 

equipped with (homogeneous) Neumann boundary conditions in and 1. Corresponding to this, 
recall the definition of the fractional Sobolev spaces H e (0, 1) as D((-A)2) or for 6 > with 
the equivalent Sobolev-Slobodeckij norm, cf. |17] . 

Define in particular the space (V, ||-||y) as H 1 with norm \\u\\ v '■= Jq |d£ii(£)| 2 d£ and denote by 
V* its topological dual, thus we have the Gelfand triplet V H V* with continuous and 
dense embeddings. As auxiliary spaces we define the product spaces V := V X H, % '■= H X H 
with V H V* continuously and densely. In the following, we define the abstract drift and 
diffusion operators mapping V to V*. 



It is well known that A can be uniquely extended to A : V — > V* via an integration by parts 

V*(Au,v) v = [ dfu(£M£) d £ = -/ H)^(^^IWII4 

Jo Jo 

The functions 4>i,4>2 ■ [0, 1] xlxl->i are continuous and satisfy some one-sided Lipschitz 
and dissipativity, as well as some (local) Lipschitz and growth condition specified below. 
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Assumption 1. There exist constants L, M, 7 > 0, /3 G M and 1 < m < 3 such that 

V?>2(6> «l,Wl) - 4>2{t,2,V2,W2)J \W1-W2J V 7 

^^'"lVf U> )</3H 2 + iH 2 -7H- +1 , (A2) 

M£>*>«0I <M(1 + H m + H), (A3) 
Mt,v,w)\ < M(l + \v\ + \w\) (A4) 



(A5) 



and moreover 

|0i(£i,ui,wi) " 01 (6,^2, ^2) | < 1(1 + M ro_1 + H 1 "'" 1 ) 

• (16 - 61 + \n - v 2\ + \wi - w 2 \). 

\<t>2(£l,Vl,Wl) ~ <fa(&,V2, w 2 )\ < L(|6 - 61 + \vi - v 2 \ + \wi - w 2 \). (A6) 
for all 6 6)6 e [0> l]i v i v\,V2, w, wi,w 2 € M. 



A typical form of 0i is a polynomial in v of odd degree and negative leading coefficient together 
with a linear perturbation in w, for example cj>i(v,w) = v — v 3 — w similar as in the FitzHugh- 
Nagumo system described in Section 6. With these assumptions we can define the Nemytskii 
operator <J> : V — > V* by 

{*(v,w)){Q .- [ h [e A e )M0) ) 
and estimate for vi, v 2 € V, Wi, w 2 € H with the embedding V t — >■ L°°(0, 1) as follows 

v *(^(v 1 ,w 1 ), (v 2 ,w 2 ))v < C(\\v 1 \\y\\v 2 \\ v + {1+ || wi\\h) \\v 2 \\h + 0- + \\m \\h + \\vi \\h) \\w2 \\h)- 

Now that the assumptions on the drift of equation ([T| are stated, let us focus on the noise process 
rj, formally given as the derivative of some Q- Wiener process y/QW(t), where (W(t))t>o is a 
cylindrical Wiener process on w.r.t. an underlying probability space (fi, J 7 , Ji,P). In order 
to obtain a solution to ([I]) (in a sense specified below) we cannot treat the case of space-time 
white noise, i.e. Q = I, but only the case of colored noise with nuclear covariance operator Q. 
The assumptions on Q and b are combined below. 

Assumption 2. Let 6 6 (h, l) and Q G L(H) be symmetric and positive definite. Assume 
that tr(— .A) e Q < 00, in particular Q admits an integral kernel of the form 

{Qv)(0 = JJ q(M)q(j>X)v(j>)dCdp, geH e ((0,l) 2 ). (A7) 

The scalar valued noise intensity b : [0, 1] x M — > R is supposed to be bounded and Lipschitz 
continuous, in particular 

|&(M|<1 and m 1 ,v 1 )-b(&,V2)\ < (16-61 + \vi-v 2 \) (A8) 

for all 6 6)6 £ [0) 1] an d f, 1*1,^2 € K. 
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In the same manner as for the nonlinear drift, we define the Nemytskii operator f£?(-u)j(£) := 
b(£,v(£,)). Furthermore, let x = (v,w) T G V and 



A(x) :-- 



Av 




+ $(v,w), B(x) :-- 



B(v)VQ 




Then, the stochastic evolution equation for X(t) := (v(t), w(t)^ T corresponding to equation ([I]) 
is given by 

dX(t) = A(X{t))dt + B(X(t))dW(t), X(0) = x . (2) 

Existence and uniqueness of a solution for such an equation can be studied in the variational 
framework, see for example |16| . in particular a recent extension in |14| . The following theorem 
is a corollary to |14| Theorem 1.1]. 

Theorem 2.1. Let T > 0, x G L 2m (n, F, P; H) and (|Al) -(|A4| , jA7| ), jA8| ) are satisfied. 
Then, equation ^ has a unique variational solution (-X"(£))tg[o,T] which satisfies 



E 



sup ||X(i)||^+ / T ||X(t)||2 d t 
te[o,T] jo 



< oo. 



Proof. We have to verify the conditions (H1)-(H4) in |14j . Of course, for X\,X2,X3 G V the map 
s i — y v*(A(xi + sx 2 ),X3)y is continuous on M, hence (HI) holds. Furthermore, we can easily 



obtain the monotonicity of A by the one-sided Lipschitz condition (Al) 



2 v *(^4(xi) -^4(x 2 ),xi — x 2 )v < — 2||«i — U2II + 2L||xi -x 2 |' 2 



Moreover, B is Lipschitz continuous in the Hilbert-Schmidt norm since 

||B(xi) -B(x 2 )\\ ns < \\B( V1 ) - B(v 2 )\\ LiH e >H) \\(-A)I^Q\\ nS 

< c e \\ Vl - v 2 \\ H {tT{-A) e QY . 

The latter inequality holds because of the embedding H s (0, 1) «— > L°°(0, 1) for > i with 
constant c#. Thus, we have verified (H2) with p = 0. ( |A2[ ) together with the fact that 
\\ b {x)\\l(h) < 1 implies 

2 v *{A(x),x) v + ||£(x)||^ 5 < -2||s||£ + 2max{L,/3 + l}||z|& + trQ 



for all x G V and (H2) holds with a = 2. With ( |A3[ ) and (A4) the matching growth condition is 

\\A(x)\\ v * < \\x\\ v + sup f \<h(ZA£)MZ)M0\<%+ ^p [ \<h(ZM0M0Ht)\<% 

\\u\\ v =lJo \\u\\ H =lJ0 

<||x|| v + 2M sup ! (l + \w{0\ + \v(0\)H0\^ + M sup f \v{£)\ m \u{0\dZ 

\\u\\n = lJo \\u\\y=lJo 

< \\x\\l + 2M(1 + \\w\\ H + + ciM||«||v ^ jf K£)| 2m_2 d£J , 

where we used the embedding V ^ L°°(0, 1). The latter is finite if m < 3, hence 

P(x)||^<C(l + ||x||2)(l + ||x||^- 2 ) 
and condition (H4) in |14] holds. □ 
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3 The Finite Differences Scheme and the Main Result 



In this section we briefly describe the approximation scheme and state the main convergence 
result in Theorem 3.1 Equation ^ is evaluated on a simple grid approximating the domain 
(0, 1), given by ^{0, ...,n} and the vectors v n ,w n G M n+1 representing the functions v and 
w on this grid. Furthermore, we introduce the spaces V n = M ra+1 with the discrete boundary 
conditions u\ — uq = u n — u n _i = corresponding to the homogeneous Neumann boundary. We 
equip V n with the semi-norm := n Y^k=i( v k ~ v k-i) 2 - In order to imitate the variational 
approach in the discrete setting we want a summation by parts formula to hold, i. e. 

n n 

h E ( Anv )k Uk = - n ^Z( v k - v k -i){u k - Uk-l) (4) 

k=0 k=l 

for all v G V n , u G M n+1 . The resulting finite difference approximation A n of A is just the 
discrete Laplacian, given by 

(A n v) k := n 2 (v k+ i - 2v k + vjfc-i), l<k<n-l, v G V n 

together with the appropriate boundary values (A n v)o = n{v\ — Uq), (A n v) n = n(v n — i> n — l) 
which are for v G V n . The nonlinearities are simply evaluated point-wise, i. e. 



<f>i !2 (v,w) := Ui,2(~,v k ,w k )\ 



0<k<n 



and b n {y) := diag (b 



,v k 



0<k<n 



for v,w G V n where by diagv we denote the diagonal matrix with the entries of the vector v on 
the main diagonal. In the next step let us construct the approximating noise in terms of the 
given realization of the driving cylindrical Wiener process (W(t))t>o- Recall that 

(5) 



(W(t), y/nhh^ fcO =: /Sfc(t), 1 < k < n 

L n 1 n J 

defines a family of n iid real valued Brownian motions. The spatial covariance structure given 
by the kernel q is discretized as follows, 



k-l I l-l 



q(^C)dCd^ l<k,l<n, 



together with q^ l := 0, 1 < I < n. With this discrete covariance matrix we can replace q by a 
piecewise constant kernel, i.e. for 1 < k < n and £ G [^p, r) 



?(£,CW*,C)dC 



n 

- E <i 
i=i 



r i-l l_- 



1 = 1 



^E<^)= : (v^wt)) fc , 



1=1 



where P n h = ((h, y/nlri-i i\))i<i< n - Denote by W n (t) = P n W(t) the resulting n-dimensional 
Brownian motion and the finite dimensional system of stochastic differential equations approxi- 
mating equation ^ or rather ([I]) is 



dv n {t) = A n v n (t) + ^(v n (t),w n (t)) dt + b n (v n (t))^/Q^dW n (t), 



dt. 



(6) 
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Standard results on SDEs imply the existence of a unique strong solution (v n (t), w n (t)} to (|6]), 
see e. g. |16[ Chapter 3]. In order to formulate the convergence result, we embed the processes 
v n and w n into the space C([0, T], V) by linear interpolation w. r. t. the space variable £. Given 
v & V n , define 

m ■= n-k+ i) Vk + (k— <H_i, e e *] • 

Denote by i n : — )• V; w i— > v the embedding given above. Furthermore define 



X n (t) :-- 



v n (t) 
w n {t) 



and E n (t) := X(t) - X n (t). 



Theorem 3.1. Suppose Assumptions 1 and 2 hold, vq S C°° (Q, J-, P; C([0, 1])) and wq £ 
J 7 , P; C 1 ([0, 1])). T/ien /or every 1 < p < oo there exists a finite constant C = 
C(vo,wo,p,T,Q,L,f3,~f) such that 



E 



sup \\E n (t)\\ P n 



< C IE 



\E n (0)\\% 



+ E 



I n {v(t)) 2 dt" P 



2 \ 2p 



je[o,T] 

which converges to as n — > oo. The rate of convergence is implicitly given by the expression 

n .k 2 i 

*»(«(*)):=(£ ^^(t/Ct.O-t/Ct.Cii)) dC) 1 . 

fe=l n 

// m addition v <E C p (fJ, J 7 , P; L 2 ([0, T], /or some a > and 1 < p* < oo, f/ien the rate 

of convergence is explicitly given by 



E 



sup \\E n (t)\\ p n 
*e[o,T] 



< C E 



|^(0)||^ 



2p + n -|min{l,a} 



/or all p < p* and some constant C = C(vo,wo,p,p* ,T,Q, L, ^,7, |H|£ P * (l 2 ([o,t] ,#!+«)))• 



4 Uniform A Priori Estimates 



In addition to the a priori estimates in Theorem 2.1 the proof of Theorem 3.1 requires such 
estimates corresponding to the approximating solutions v n and w n . These estimates are formally 
obtained by applying Ito's formula to the square of the "H-norm and we will do the same in the 
discrete case and obtain an exponential a priori estimate for the L 2 ([0, T], y)-norm of v n and 



L°°([0, T], V)-norm of w n uniform in the parameter n. Moreover due to (A2), we have control 
over the discrete l m +i-norm of v n uniformly in n. For this purpose define 



00. 



fc=i 
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Lemma 4.1. Let n E N. The approximating solution (y n {t), u>"(t)) satisfies the estimate 
E 



q{v n (t))+q(w n (t))+2 \\v n {t)\\l&t + 2 1 / c + iK(t))dt 

Jo J (7) 

< e 2max{L ^ }T (E[ZJ(u n (0)) +/2(w"(0))] +TtrQ) /or aHte [0,T]. 

Proof. Ito's formula applied to (f/?) 2 + (^fc) 2 , the summation by parts formula Q and the 
dissipativity condition on <pi,<p2 hi (|A2|) imply 



5K(t))+iJ(« n (t))+2 /V«ll»d* + 27 fc +1 {v n {s))ds 

Jo Jo 

<q(v n {0)) +q(w n (0)) +2L [ q(v n {s))ds + 2p [ q(w n (s))ds 



'0 JO 



+ 

n 



- 70 fe,/=i v - 70 fc,/=i 



With a standard localization argument the stochastic integral vanishes after taking the expecta- 
tion and since b is bounded we obtain 



E 



t rt 

ni „\ [■ 2 



q(v n (t))+q(w n (t))+2 \\v n (s)\\ 2 n ds + 2 1 C +1 (v n (s))ds 

Jo Jo 



< E 



%(v n (0)) +l%(w n (0)) +2max.{L <j> ,P} / q{v n {s)) + q{w n (s)) ds + itr Q 







Gronwall's lemma now yields the result. □ 

Unlike u, the equation for w has no regularizing linear part. Nevertheless, one can improve the 
a priori estimate significantly if the initial condition has more regularity, because the coupling 
with v is Lipschitz continuous. 

Lemma 4.2. The following improved a priori estimate holds for every n E N 



E 



sup \\w 

te[o,T] 



<£ e 5^«{£>/3W E [|| w n( )||2 + £(„»(())) +%{w n (0))] + TixQ\ (8) 



Proof Consider the equation for w^(t) for itself and a straightforward calculation yields 

\wm-™%-i(t)\ 2 = \™m-™um 2 

+ 2^K( S )-<_ 1 ( S ))(0 2 (|,^( S ),<( S ))-0 2 (^,^ 1 ( S ),<_ 1 ( S )))d S 







<3L -<_ 1 (,)| 2 d, + 2L A^ + KW-^WI 2 )^. 

jo jo 



S 



Hence, a summation over k and an application of Gronwall's lemma yield 

i-T 



E 



sup \\w n {t)\\l 



<e 3LT \\w n 



+ 2LE 



te[o,T] 

Now use the estimate ([7]) to conclude ^ 



l + \\v n 



dt 



o 



□ 



By a similar strategy as in Lemmas 4.1 and 4.2 we can also obtain exponential a priori estimates. 



This is due to a standard decomposition of the exponential of a martingale into a supermartingale 
and the exponential of its quadratic variation process. The latter can be bounded independent 
of the size of trQ by the estimate in l m +i- The precise statement is the following. 

Proposition 4.3. Let v G C°° (Q, T, P; C([0, 1])) , w G £°°(Q, J^P; C^flti, 1])) and a > 
arbitrary. Then 



E 



exp|a( / (||-5 n (i)||^ + 7C + l(^)))dt+ sup ||^(t)|| 2 y 
o *e[o,T] 



< 



uniformly in n G N, where C fe^i = C (a, vq, wq,T, Q, L, to) is an explicitly known finite constant. 



Proof. Note, that the ltd formula in the proof of Lemma 4.1 and Gronwall's lemma imply 
T {\\v n (t)\\l + 7C+iKW)) dt < ie 2max ^> T fe(^(0)) + £( W »(0))+TtrQ + M T 



o 



where My denotes the stochastic integral. Also, we have seen in Lemma 4.2 that 

r-T 



sup \\W"{i < e 3LT ( \\w n 



+ 2L/ \\v n (t)\\ 2 n dt). 



te[o,T] 

Hence, a combination of the two inequalities above yields that for any a > 

exp(a( f ' {\\v n (t)\\l+7C+l{v n (t)))dt+ sup \\w n (t)f n 
\ x Jo te[o,T] 

< exp |aCi(^(/(0))+^(w n (0)) + |K(0)|| 2 + Ttr • exp (aC 2 M^j. 



O) 



Here, d = (2L + l) e 5max{L,/3}T and C2 = e 2max{L,/3}T + 2 Le 3LT+2 ™x{L,/3}T Now recall that 

for every continuous local martingale vanishing at t = the process 



a 



exp[aM t -^-(M) t 



is again a continuous local martingale, thus a supermartingale by Fatou's lemma and therefore 
E[Zf ] < 1 for all t. With this information we can derive 



E 



exp ( -M t 



< E 



(10) 
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Now we can take the expectation in ([£]), Holder's inequality and then (10) implies 
E 



exp|al / {\\v n {t)\\ 2 n + 1 C +1 (v n (t)))dt+ sup\\w 
o te[o,T] 



< E 



< E 



exp ^2aCi(^(v n (0)) +q(w n (0)) + \\w n (0)\\ 2 n + T tr Q 
exp [2aC 1 (q(v n (0))+q(w n (0)) + \\w n (0)\\ 2 n + TtrQ 



E 


exp | 


E 


exp | 



.,2^2/ 



The quadratic variation (M)t can be bounded by the integrals on the left hand side, in particular 
by the Z^ 1+1 -norm of v n . 



^ = 4 f E (Ks.«?(*))^(t)^) 9 d* 
j ° k,i=i 



n 3 .lo 



E K(*)9fc,z) d* 

c,/=l 

. .-j n n 

s^/ (E(« ! )(E(*)> 

1/0 fc=l k,l=l 

<4trQ [ T q(v n (t))dt 
Jo 

pT m + l 

<4trQ ^ 5C+1 (^(t)) +5 «^if^_2_ T N )- 1 ( !/, 



2 \ 8 (m+l) J 

where we used Young's inequality in the last line for arbitrary 5 > 0. In particular, we can 
choose 5 = 7(32aCf tr Q)~ l such that 



E 



exp(a( f T (\\v n (t)f n + 7 l^ +l (v n (t)))dt+ sup \\w n (t)\\l 
\ V J0 te[o,T] 



< e IC 6 TtrQ E 



exp 2aC 1 (iJ(t; n (0)) + #(u; n (0)) + \\w n (0)\\ 2 n + Tti Q 



(11) 



where is the constant appearing in Young's inequality above. Finally, observe that for v G V n 
the piecewise linear interpolation v is in V and its weak derivative is explicitly given by the 
differences n(vk — Vk-i)- Hence 



l*l| 2 v 



k=l n 



if &£ = ^E { v k ~ v k -i) 2 - 



k=l 



Moreover, note that for vq G C([0, 1]) and wo G C 1 ([0, 1]) we can define the point- wise evaluation 
v n (0) and w n (0), which satisfy 



£(w*(0)) < lbo|| 2 c ([ o,i D and / 2 n K(0)) + [|uff(0)||f, < [|wo||^ ([0>1 ]) 



and we found an upper bound for (11 ) uniformly in n G N. 



□ 
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5 Error Estimation, Proof of Theorem 



3.1 



This section is devoted to the proof of Theorem 3.1 and the estimation of the error in 
£ p (p, T,F;C([0,T],H)). We proceed in several steps, first some path-wise estimates and 
then using the Burkholder-Davis-Gundy ineq ualit y together with the a priori information on 
the approximating solution X n in Proposition 4.3 to prove that X n converges strongly to the 
exact solution X. 



Step 1: Ito's formula applied to the square of the %-norm (which is available for the variational 
solution) implies 



\E n 



\E n 



v* (Av(s) - L n A n v n (s),v(s) - v n {s)) v ds 

+ 2 / v *(<S>{v(s),w(s)) - L n <l> n {v n (s),w n (s)),E n (s)) v ds 
Jo 

+ 2 J (v(s) - v n (s), (b(v(s)) ,/Q - t n b n [v n (s)) y^Pn) dW(s)) 

+ f \\B(v (s)) - t n b n {v n (s)) y^Pnfm ds 
Jo 



E\ + Eo + E3 + E4 



(12) 



Each of the error terms Ei,E2, E4 can be split into two parts with the first one only involving the 
difference E n and the second one the approximation of the parameters of the equation. For the 
stochastic integral in E3 we use the same argument after applying the Burkholder-Davis-Gundy 
inequality later. Now let us begin with, for instance E\ where an integration by parts yields 

Ex = -2 [ t \\ v ( s )-v n (s)\\ 2 v ds + 2 f v *(Av n {s) - L n A n v n (s),v{s) - v n (s)) v ds. 



Moreover the one-sided Lipschitz condition (Al) for the nonlinear drift part and the Lipschitz 
continuity of B in ([3| yield 

E 2 < 2L f t \\E n (s)\\ 2 n ds + 2 f v *{${v n (s),w n (sj) - L n $ n {v n (s), w n (s)) , E n (s)) v ds 
Jo Jo 

E 3 <2c e tx(-A) e Q [ \\E n (s)\\ 2 n ds + 2 [ \\B(v n (s)) - i n b n (v n (s)) v^nllls <±s. 
Jo Jo 

In the following, we estimate each of the error terms for fixed s, thus we drop the time dependence 
in the notation. 



Step 2: Approximation error of the Laplacian 



At first, we will prove that the error term coming from the linear part converges to if v G V, 
i. e. we need only the guaranteed regularity of the variational solution. The rate of convergence 
is not uniform and given implicitly in terms of the solution v. Secondly, if we assume additional 
regularity on v, we deduce an explicit rate in terms of n. Let <p EV, then 

v , (Av n - i n A n v n , <p) v = v * (Av n , <p) v - (t n A n v n , if) 
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since the linear interpolation is in V. In the first term an integration by parts yields 

n 

v*(Av n ,<p) v = -» J>J - - (13) 

k=l 

In the second term we use the summation by parts formula Q, hence 

fi h. h. 

{i n A n v n , <p) = J2 (A n v n ) k (n£-k + 1M0 d£ + (A n v n ) k _ 1 (k - d£ 

k=l n n 

n 

= £ + Vi^-i 

k=l 

n 

= -n 2 Y,K-Vk-i){(9k-g k -i) + (9k-9k-i)), (14) 

k=l 

where we denote the integrals by g^ and gk-i, respectively. Note that the boundary term vanish 
because of the boundary conditions for v n . The next step is to replace <p by E n but since 
the computations differ for v and the piecewise linear v n we split these terms and obtain the 
following lemmas. 

Lemma 5.1. For all n £ N it holds that 

V*(Av n -i n A n v n ,v n ) v < 0. 

Proof. Since v n is piecewise linear we can compute g^ and g^, which are 
9k = ±vS + ±<V-l and Sk = l^ 1 + 1^. 



Obviously, we also have v n (^) = and ^ n (^p) = v k-i- Equations (13) and (14) now read as 



n _ n 

v *(Av n - t n A n v n ,v n ) v = -n^K - ^„i) 2 + y J>J? - vt-i? 

fc=i k=i 

n n— 1 

fc=2 fc=l 

With Cauchy's inequality we can bound the latter terms by the former ones, thus the statement 
is proven. □ 

Lemma 5.2. Let tp € V , then 

v *{Av n - L n A n v n ,<p) v < 2\\v n \\ v ■ I n (<p), 

where 



r 2 - 

fc=l n 
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Proof. The factor with tp in equations ( 13 ) and ( 14 ) can be written as follows 
I - - (9k ~ 9k-l) + (9k ~ 9k-i) 

k 

Zi ^~ k + ^ [>(*) - " (v(*ir) " viS ~ 

+ (k- nO [[ V (0 ~ <p(*=±)) ~ + I) - ¥>(*))] d£ == Cn(<p), 

because (n£ — /c + 1) + [k — n£) = 1 and the integrals appearing in g\--\ and g^ are shifted to 
the same interval. Rewrite c n (tp) with the fundamental theorem of calculus and use Jensen's 
inequality to obtain 

(k k 
J^H-k+i) | n ^(o-^(c-i)dcde 
n ^ 

v 9 

e 

fe-i 



+ 



< 2rT 



7 



n / 

yccW(c±£))V- 

All there remains to do is an application of the Cauchy-Schwarz inequality. 



(15) 
□ 



Lemma 5.3. Let I n (tp) be defined as in Lemma 5.2. 

i. If tp £ V, then I n (tp) < 4||</j||y and moreover lim ri _ >0O I n (}P) = 0. 

ii. If tp 6 H 1+a (0, 1), the explicit rate of convergence is given by I n (<p) < 2 a+2 n -a ||(/j||#i+a . 

Proof. At first, we will consider the second assertion. Thus, let tp G C 1+a ([0, 1]) and denote by 
(p' its first derivative. Recall the definition of c n {tp) in the proof above. An additional integral 
in £ over the same interval yields another n, while we can insert </>'(£) m t° the square above. 
The distance \Q ± ~ — £| is of course bounded by ~, hence 



CnM 2 < 2 



2a+3 -3-2q 



ti / n 



WWCfl) 2 , (^(fl-^(C±£)) 2 

,. , . „ , ic-ei 2a+l 1 \c±^-e a+1 



and we have shown ii., since the sum of the latter expression is uniformly bounded by the 
Sobolev-Slobodeckij semi-norm of tp in 

In a second step, we prove that if tp E V, the expression I n (tp) is bounded uniformly in n 
and we can indeed approximate tp by sufficiently smooth functions in V to obtain the desired 



convergence. For this purpose, we estimate ( 15 ) simply by 



c n (<p) 2 < 8n 



-3 



l^(C)l 2 + l^(C±^)l 2 dC, 
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hence J n (y) < 4[|^[|v. Now, let e > be given. Clearly, C 1+Q ([0, 1]) C V densely for a > 0, 
therefore one can find (p £ G C 1+a ([0, 1]) such that \\ip — ip e \\v < §• Furthermore, we can find 
no 6 N such that for all n > hq 



If we turn our focus back on the original error term E\ and apply the previous lemmas, we have 
shown that ^ ^ 

E x <-2( \\v(s) -v n {s)\\ 2 v ds + 2 f \\i n (s)\\ v I n (v(s))ds. (16) 
Jo JO 

Step 3: Approximation error of the nonlinearities 



The nonlinear terms in E2 are estimated with the (local) Lipschitz conditions (A5) and (A6) 
together with the uniform a priori estimates on the approximating solutions v n and w n . The 
error term in E2 is 

y. «>"( S )) - L n <S> n (v n (s),W n (s)),E n (s)) V 

which consists of contributing parts of fa and fa. Since the assumptions on fa are much 
stronger, we exemplarily estimate the one for fa in the following. 



k=l 



~{ni-k + I) fa (£, - (ft - nC)0i(fci, i£_ l9 " v n (0) ^ 

^E/Jl [iK-fc+i)i-i^(c,^(e),* n (0)-^i(i<,Oi 

A;=l n 

+ |ft - <| • |fc (£, 5»(0, tS n (0) " & l] Kf ) " v n (0\ ^ 

fc = l n 

■ [i + k - vi_,\ + \wt - wt_i\j\ ho - v n (o\ 



where we used (A3), |£ ra (£)| < \v%\ + \v? | and that \v n (^) — = (k — n£)\v% — as 



well as |v n (£) - v%_ x \ = « - fe + - v n 



fc-ll- 



Of course, the same holds for w n . With 



Cauchy-Schwarz' and Young's inequality we can further estimate 

n n ,.k 

<^E( 1+ Kr -1 ) 2 + f^E (i + k - ^ n -ii + k - <-iD 2 L>(£) - ^)i 2 d ^ 

k=l k=l n 

<±l(i + C+i(0) + + + W^Wl) ll« " v n \\%, 

since 2(m — 1) < m + 1 for m < 3 and the gradient of the linear interpolation is given by these 



finite differences, see Proposition 4.3 For simplicity we estimate fa's part of the error in the 
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same way, hence 

E 2 <6L r(l+||^( S )||2 y + ||^( s) ||2 y )|| jB n (s) ||2 {ds+ i L f\l +C+1 ( v n {s) )) ds . (17) 

Jo Jo 
Step 4: Approximation of the covariance operator 

The Ito correction term gives us a Hilbert-Schmidt norm which, in our case, is given by the 
L 2 -norm of the associated kernel. Thus, we can write this part of the error term E4 similar to 
the one of the nonlinear part in Step 3. Recall that P n h = {{h, \fn\,i-i i_\))\<i< n - 

\\B(y n )^Q-t n b n (y n )^P n \\l s 

t n 

' L 77, 1 77, J 



V 

Jo 



k=i 



« - + vl) ql^^iC) + (k — ^-1) £ ^_i^[^,i)(C)l dCd£. 

* L 5 77, I * L 77, ' 77, V J / 



A closer look reveals that the error consists of parts Si contributed by the approximation of b 
and S2 by the approximation of the kernel q. Therefore we split the norm square into two parts 
where the first one is 

2 E / ^{K^^-H-k + l^v^-ik-nOb^^)) q^CfdCdt 

k=l J ° J ^rT 

Similar to Step 3 for the nonlinear drift term, since b is Lipschitz continuous, 

(b{t, v n (0) ~H-k + l)b( k -,v n k ) -(k- nO&(*=l, ^_i)) 2 <U& + K " 
and the summation over n yields 

Si < £trQ(l + \\v n f v ). 
The second error 5*2 term due to the approximation of the covariance kernel is given by 

n h_ _L 



fc,/=l 



n n 



where we already used (n£ — fc + 1) < 1 as well as (k — n£) < 1. Exemplarily we do the estimate 
for the part with q^ l but note that the procedure remains the same in the one with q%_ 1 ; except 
for a constant. Since b is bounded, we obtain 

n k I 

?1 _ /c / k I 

= 4 E c) - - 2 /: q {p, v) ^ ^ 2 d C d^ 

k,l = l n n n n 

n k 1 fc 1 

* 4 - 2 £ /I /I £ /: c) - .)i 2 de dc dp d„. 



fc,J=l • 



n n n n 
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Since |£ — p\ < - as well as \C t — 77 1 < - we can expand this expression to the Sobolev-Slobodeckij 
semi-norm of q. 

<2^n-™ [fff 1 Igfr0-g(p,q)| 2 dgdCdpd, 

yyiio de-H 2 +ic-r?i 2 ) e+i 

<2 e+3 n- 2 lg||^. 

In summary, we have shown that (the modified constants are due to the omitted terms) 

E4<2c e tr(-A) e Q I \\E n {s)\\^ ds + | tr Q [ (l + \\v n (s)f v ) ds + 2 3e+5 n- 29 \\q\\ 2 H e ■ t (18) 
Jo Jo 

Step 5: Control over the supremum 

The four previous steps allow a first estimate non-uniform in t. For this purpose define the 
processes 

L n (t) := 6L(1 + ll^^ll 2 / + \\w n (t)f v ) + 2c e tr(-A) 9 Q 

and 

K n {t) := i(4L(l + C + iKW)) + 2trg(l + [|w"(<)||^) 

+ 2 39+5 ||g||^)+2||r(t)||yl4t;(t)). 



Now we apply Ito's formula with the function f(x) = x p , p > 1 for the square of the %-norm of 
the error E n (t). By Steps 1-4 

\\E n (t)\\%<\\E«(0)\\%+p [ t L n (s)\\E n (s)\\^ds+p f K n {s)\\E n {s)\\%- 2 ds 

Jo Jo 

+ p [\\E n (s)\\%- 2 dM(s) + p -^ f* \\E n (s)\\ 2 ^~ 4 d(M) s , 
Jo Jo 

where the local martingale M(t) is given by 

M(t) := 2 J\v(s) - v n (s), (B(v(aj)y/Q- t n b n {v n (s)) ^ P n ) dW(s)). 
Its quadratic variation is simply 



(M) t = / \\v(s)- v n (s)\\ 2 H \\B(v(s))^Q- L n b n [v n (s))^P n \\ 2 m ds. 
Jo 

Obviously, there appears exactly the same term as in E4 which we already estimated point-wise 



in t. Thus, use (18) and Ito's product rule to conclude 

sup e-^ L p( s ) ds \\E n (t)\\ 2 ^ < \\E n {0)\\% + / e-J?W d '/C^(t)||^(*)||g~ 2 dt 

,-t 



tew jo (ig) 



+ sup f e~ to ' L > )dT \\E n {s)\\ 2 £- 2 dM{s). 

t£[0,T]J0 
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The processes L™(t) and Kp(t) are given by 

L n p (t) := P L n (t) + p(p - l)c0 tr(-A) e Q 

and 



K;{t) := pK n {t) + Ip(p - 1)( tr Q(l + IK(i)l^) + 2 3 *+ 4 ||< Z " 2 



The first integral in ( 19 ) can be estimated with Young's inequality as follows 

[ T e -ti L ^ ds K"(t) \\E n (t)\\%- 2 dt< sup e-ti L ^ ds \\E n (t)\\%~ 2 C ' K?{t)di 
Jo " te[o,T] Jo 



<ti sup e-tiW^E"(t)\\%+l-( f K; 
p te[o,T] p \Jo 



(t) dt 



Now, take the expectation in (19) and the Burkholder-Davis-Gundy inequality bounds the 
supremum of the stochastic integral from above by its quadratic variation, more precisely 



E 



sup e -/oW ds ||i^(t)g 
*e[o,T] 



< pE 
+ 3pE 



E' 



n (C\\\\ 2 P 

H 



+ E 



T \ p ' 

n 
P 



K$(t) dt 



- 2 /oWds[|£» (i) [|4 P -4 d(M)i 



(20) 



The quadratic variation is again estimated as above in terms of the left hand side and K™ (t) . 
This is done with Young's inequality as follows 



3pE 



T 



< E 



2p-l ft 



So L p( s ) ds \\ E™|V\ll 2p_1 



sup e 2p Jo-P^-\\ E r. yuJm 
te[o,T] 



< 3(6p-3) 2p - 1 E 



■ y £ ti L ^ ds \\B(v(t)) - Lnb n {v n (t)) VO^Pnfus dt 
J^e-Uo^^llBivit))^ - i n b n {v n (t))^PJl s dt) 



+ §E 



sup e-ti Ln ^ As \\E n {t)\\ 2 £ 
ie[o,T] 



The latter is absorbed by the left hand side and the former is again essentially E4 which we 
already estimated in Step 4. With (18) and (20) follows 

E 



sup e-ti L » ds \\E nf ^ 21 ' 



*S[0,T] 



H 



< 2pE 
+ C 4 E 



E n (0)\\% 



+ C 3 E 



T 







K£(t) dt 



sup e-^ L ?W dr ||^ n (s)||Jdi 
se[o,t] 



17 



where C 3 := 6(6p - 3) 2 p _1 2P + 2 and C 4 := 12(6p - 3)* -1 2Pcb tr(-y4) 9 Q. Now, GronwaU's 
inequality implies for all p > 1 



E 



sup e~fo L p( s ) ds \\E nu ' 112 '' 



H 



< e 



C A T, 



2p& 



|E"(0)[|* 



+ C 3 E 



K n (t) dt 



• (21) 



_ie[o,r] 

This preliminary error estimate yields the desired one via Holder's inequality, since we have 
(21 ) for every p > 1 and can control the exponential by Proposition 4.3 Therefore, we have for 
p > 1 



E 



sup \\E n (t)\\ p H 
te[o,T] 



< E 



< E 



sup e -2/o^Wd S || jB n (t) ||P e |/ ( ri«( t ) 



*e[o,T] 



sup e-Zo^W^H^^)!!* 



te[o,T] 

<C^e c * T \2p& \\E"{V)\\% 



E 



J T L"(t) dt 



^p 



K$(t) dt 



j_ 



(22) 



It remains to study the convergence of K n to 0. For this purpose, we fix p = 1 at first. Then 



iff (f) = K n {t) and by Proposition |4~3 

fT 



E 



K n (t)dt 



< Cn' 1 + 2E 



\v n (t)\\ v I n (v(t))dt 



Since by Lemma 5.2 I n (v(t)) < 4||u(f)[|y and both v and u n are (uniformly) bounded in 
C? (p,, J 7 , P; L 2 ([0, T], V)) we can use Lebesgue's dominated convergence theorem to deduce 

fT 



E 



K n {t) dt 



0, hence E 



sup \\E n (t)\\ H 

te[o,T] 



as n — > 00 



by Lemma 5.2 Thus, the exponential moment estimates carry over to the limit v and we can 
apply the dominated convergence theorem in cases p > 1 which concludes the first assertion of 
Theorem 3.1 If moreover v G C p * (Q, J 7 , P; L 2 ([0, T], H 1+a )"j for some a > and p* > 1, then 
again by Lemma 5.2 

cT 



E 



K n {t) dt 







< Cn~ p + C'E 



< Cn" p + C'E 



< Cn~ p + C"n- ap E 



\v n (t)fydt 



In{v(t)) 2 dt 



\v n (t)\\ldt^ 2(p *- p) 



dt 



p* 



E 



I n (v(t)) 2 dt 



V[t\\ H U 



hence the second assertion is proven for all p < p* . 



□ 
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6 Applications 



As an application for our results in Theorem 3.1 we consider the spatially extended FitzHugh- 
Nagumo system with noise. Originally, this was stated by FitzHugh [3j as a system of ODEs 
simplifying the famous Hodgkin-Huxley model [8] for the generation of an action potential in 
a neuron in terms of a voltage variable v and a so-called recovery variable w. Its spatially 
extended version is a model for the propagation of the action potential in the axon of a neuron. 
See e. g. the monographs (2j [12] for more details on the deterministic case. Now consider this 
system subject to external noise only in the voltage variable v. Together with the original 
parameters from [J) this reads as 



v(t, = a%v{t, + v(t, - -v(t, 3 - w(t, + r)(t, o, 
w(t, f ) = 0.08(u(t, £) - 0.8w(t, f ) + 0.7) , t > 0, f e (0, L) 



(23) 



equipped with homogeneous Neumann boundary conditions in and L. One can immediately 
see, that (23) is of the type Q. The first mathematical rigorous analysis of this equation in 
the context of mild solutions can be found in [Tj. It has been observed in e. g. |18j that this 
system has traveling pulse solutions (Figure [T]) which may break down due to the influence of 
the external noise, hence there is no transmission of the signal from to L. This phenomenon is 
usually referred to as the propagation failure and one is interested in calculating its probability 
depending on the strength of the external noise, see [18] for a heuristic approach. 



I a: 




100 120 140 160 



ee[o,i] 



Figure 1: A traveling pulse solution of (231 obtained 



by the finite difference scheme propagating along (0, L) 
at three different times (solid, dashed, dotted). 

An appropriate estimator detecting the event of a propagation failure is given by the integral 

$(«) := f v (£) - v* df , 
Jo 

which significantly differs in cases with or without the traveling pulse based on the observation 
from Figure [T] Here, v* ~ —1.1994 is the voltage component of the unique real equilibrium 



point for the drift in (23). Given <£•, the event of a propagation failure can be defined by 

$(«(*)) < K for some T <t<T, 
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for some appropriate threshold k > and some initializing time To. The quantity of interest is 
the probability 



Q 



min $>(v(t)) < k 

T <t<T y ' 



PQ,K 



of propagation failure depending on the noise covariance Q. Essentially we have to estimate the 
parameter of a Bernoulli distributed random variable. The sample average 



1 rn 
m — * i 



PQM V ) '— rn Z^ i {mm To < t < T $(i,('=) (*))<«} 
k=l 

based on m iid copies v := (v^)i<k<m is a natural estimator. Since the solution v is not 
given explicitly, we can approximate pQ tK based on independent numerical observations v n '■= 
({? n 'W) 1 <^< m , as independent realizations of v n only. Besides the statistical Monte-Carlo error, 
there appears additional uncertainty due to the approximation of the exact solution. Theorem 



3.1 now allows to quantify this. 

Corollary 6.1. Let pQ K := pQ )K {v n ) and e > be a priori given. Furthermore, assume that 
the solution v £ £ p (f2, J 7 , P; L 2 ([0, T], H 2 (0, 1))) for some p > 2. Then, a confidence interval for 
the estimation Pq K _ £ of pQ tK with confidence level a is given by [pg K _ £ — 7,j5q K _ £ + 7] where 

7 = ( am )^ (l + 4e~ 2 C , n~ 1 ) 2 . 

Remark. The additional regularity of the solution may be obtained in the context of mild 
solutions as in [1] if Q is sufficiently regular, since the heat semigroup in the equation for v is 
analytic and therefore maps H to D(A) = H 2 (0, 1). However, this needs further investigation. 



Proof. We can estimate the probability of Pq, k being outside of an interval of size 7 around the 
estimator pg K _ e by Chebychev's inequality with 

PQ [\PQ,k ~ PQ,k-s\ > 7] = PQ [\PQ,k ~ PQ,k + PQ,k ~ Pq,k-s\ > l] 

< PQ [\PQ, k ~ PQA > i] + [\PQ,k ~ Pq,ks\ > I] 
1 



< 



7 z m 



{sup te[0iT] [$( U ('=)(t))-$(i5".(*=)(t))]> e } > m l 



1 4 

< — + ^^E 
7 z m ^f z £ z m 



sup 

te[o,T] 



<S>(v {k \t)) - <S>(v n + k \t)) 



The additional uncertainty due to the approximation of the exact solution is the mean squared 
error of &(v n (t)). 



E 



sup ($(«(*)) - $(v n {t) 
te[o,T] v 



= E 


sup ( / 




te[o,T] K Jo 







< LE 



sup \\ v (t)-v n (t)\\ 2 H 
ie[o,T] 



< LE 



sup \\E n (t)f n 
te[o,T] 



Theorem 3.1 now implies the convergence rate of n 1 and we set 

4L 



1 



a 



+ -^^Cn- 1 . 



7 2 m 7 2 e 2 m 



□ 
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